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ABSTRACT 

Meteorological and air-quality model simulations are analyzed alongside observations to investigate the 
role of the Chesapeake Bay breeze on surface air quality, pollutant transport, and boundary layer venting. A 
case study was conducted to understand why a particular day was the only one during an 11-day ship-based 
field campaign on which surface ozone was not elevated in concentration over the Chesapeake Bay relative to 
the closest upwind site and why high ozone concentrations were observed aloft by in situ aircraft observations. 

Results show that southerly winds during the overnight and early-morning hours prevented the advection of 
air pollutants from the Washington, D.C., and Baltimore, Maryland, metropolitan areas over the surface 
waters of the bay. A strong and prolonged bay breeze developed during the late morning and early afternoon 
along the western coastline of the bay. The strength and duration of the bay breeze allowed pollutants to 
converge, resulting in high concentrations locally near the bay-breeze front within the Baltimore metropol- 
itan area, where they were then lofted to the top of the planetary boundary layer (PBL). Near the top of the 
PBL, these pollutants were horizontally advected to a region with lower PBL heights, resulting in pollution 
transport out of the boundary layer and into the free troposphere. This elevated layer of air pollution aloft was 
transported downwind into New England by early the following morning where it likely mixed down to the 
surface, affecting air quality as the boundary layer grew. 


1. Introduction 

Air pollution affects human health (Mudway and 
Kelly 2000; Dockery et al. 1993; Samet et al. 2000), crop 
yields (Booker et al. 2009; Fishman et al. 2010; Sanders 
et al. 1992), aquatic ecosystems (Moffat 1998; Morgan 
and Owens 2001; Galloway et al. 2003; National 
Research Council 1995), and climate (Ramanathan and 
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Feng 2009; Shine 2000; Fishman et al. 1979; Hansen et al. 
1997). The buildup of tropospheric trace gases and 
aerosols is dependent on air-pollution emissions, mete- 
orological conditions, and the chemical processes oc- 
curring in the atmosphere. The complex nature of 
air-pollution formation and buildup can be simulated 
with a combination of meteorological and air-quality 
models, and observations are critical in evaluating the 
models. Meteorological and air-quality model simula- 
tions provide a full picture of the conditions (if accurate), 
whereas observations are more spotty. Therefore, the 
use of models alongside observations is useful for in- 
vestigating how pollutants evolve in the atmosphere, 
forecasting air quality and the climatic impacts of pol- 
lutants, and helping to formulate air-pollution and 
climatic-change mitigation plans. 
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Previous studies have shown that sea-, bay-, or lake- 
breeze circulations can contribute to poor air quality. 
For example, along the Portuguese west coast surface 
ozone levels typically become higher when a sea breeze 
is present (Evtyugina et al. 2006). In Southern Cal- 
ifornia, the peak surface ozone concentrations on 
high-ozone days have been found to occur at the 
farthest-inland location of a sea breeze’s convergence 
zone (Boucouvala and Bornstein 2003). In Houston, 
Texas, high surface ozone episodes typically begin when 
the synoptic-scale winds transport pollutants from land 
to sea prior to the onset of a bay breeze (Banta et al. 
2005; Darby 2005). As the bay breeze begins to develop, 
stagnant conditions ensue as the winds over the water 
begin to reverse direction, allowing ozone and ozone 
precursors to accumulate before being advected on- 
shore as the bay breeze increases in intensity later in the 
afternoon. In addition, a Southeast Asian study found 
that sea-breeze circulations affect the distribution of 
smoke (Wang et al. 2013). A study during July of 2011 
showed that surface ozone concentrations are usually 
higher over the Chesapeake Bay than over upwind areas 
as a result of lower deposition rates, ship emissions 
trapped in a shallow planetary boundary layer, higher 
photolysis rates, and decreased boundary layer venting 
due to fewer clouds being present in comparison with 
the adjacent land (Goldberg et al. 2014). In the Wash- 
ington, D.C., and Baltimore, Maryland, metropolitan 
areas, a bay-breeze case study (Loughner et al. 2011) 
was analyzed in which synoptic-scale winds were west- 
erly, and it was found that 1) prior to the development of 
the bay breeze, westerly winds transport pollutants from 
urban areas out over the surface waters of the Ches- 
apeake Bay; 2) as the bay breeze begins to form, stag- 
nation develops over the bay, allowing pollutants to 
accumulate as the winds begin to change to a southerly 
direction; and 3) once the bay breeze forms, southerly 
winds over the bay transport the high concentrations of 
surface pollutants that accumulated over the bay 
northward across the coastline. The bay breeze typically 
enhances air-pollution events at Edgewood, Maryland, 
which is on the northern coastline of the Chesapeake 
Bay, making it the most polluted site in Maryland for 
ozone. In addition, it was found that, once the Ches- 
apeake Bay breeze forms, surface pollutants over land 
cannot be transported near the surface across the 
coastline but instead are transported to the bay-breeze 
convergence zone where they are lofted and then 
transported downwind (Loughner et al. 2011). 

Combining model simulations with detailed observa- 
tions at the land-ocean interface is essential for testing 
the theoretical understanding of the atmosphere as ex- 
pressed in state-of-the-art photochemical models and 


for improving knowledge of the impacts of coastal pro- 
cesses on local pollution episodes. Here, we examine 
the role of the Chesapeake Bay breeze on air quality 
and pollutant transport from the planetary boundary 
layer (PBL) to the free troposphere on 11 July 2011, 
which was an observation day during two simulta- 
neous field campaigns in the region: the Deriving In- 
formation on Surface Conditions from Column and 
Vertically Resolved Observations Relevant to Air 
Quality project (DISCOVER-AQ; http://discover-aq. 
larc.nasa.gov) and the Geostationary Coastal and Air 
Pollution Events-Chesapeake Bay Oceanographic Cam- 
paign with DISCOVER-AQ (GeoCAPE-CBODAQ; 
Tzortziou et al. 2014; http://neptune.gsfc.nasa.gov/osb/ 
index. php?section=250). 

The DISCOVER-AQ project is a National Aero- 
nautics and Space Administration (NASA) Earth Ven- 
ture mission comprising four deployments of aircraft- and 
ground-based measurement networks. These measure- 
ment campaigns are designed to improve the interpre- 
tation of satellite observations to diagnose near-surface 
conditions relating to air quality. The first of four de- 
ployments took place over the Washington and Balti- 
more metropolitan area and the Chesapeake Bay 
throughout July of 2011. Two aircraft were used, the 
NASA P-3B for in situ sampling in the lowest 3 km of the 
atmosphere and the NASA UC-12 flying at approxi- 
mately 9 km with remote sensing instruments for trace 
gases and aerosols. In addition, extensive ground ob- 
servations were used to measure air pollution at the 
surface using in situ observations and aloft using bal- 
loons and remote sensing instruments. Field observa- 
tions at Edgewood during the DISCOVER-AQ field 
campaign revealed that the bay-breeze circulation was 
evident on nine days at this location during July 2011 
(Stauffer et al. 2014). The GeoCAPE-CBODAQ field 
campaign complemented the DISCOVER-AQ by mak- 
ing ship-based measurements of air and water quality 
from 11 to 20 July 2011. 

In this study, we use results from meteorological and 
air-quality model simulations, airborne and shipborne 
in situ air-quality observations, and routine air-quality 
and meteorological observations from ground-based 
monitoring stations to investigate how the Chesapeake 
Bay breeze affects air quality. We focus on the transport 
processes that took place on 11 July for two reasons: 
1) this was the only day on which ozone concentrations 
were not increased above baseline over the bay relative 
to the closest upwind ground-based monitoring stations 
(Goldberg et al. 2014), contributing to lower observed 
ozone concentrations at Edgewood than at nearby 
monitoring sites to the west, and 2) increased levels of 
air pollution were observed aloft by in situ observations 
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on board the NASA P-3B aircraft relative to surface 
concentrations (to be shown and discussed further be- 
low). Past research on the relationship of the Ches- 
apeake Bay breeze and air quality showed increased 
pollution levels over the bay waters that were due to 
stagnation that develops over the water as the winds 
begin to change direction as a result of bay-breeze de- 
velopment (Loughner et al. 2011). In addition, lower 
deposition rates and PBL heights over the water con- 
tribute to high air-pollution concentrations over coastal 
urban waters (Angevine et al. 2004; Martins et al. 2012). 
We focus on this event to gain an understanding of why 
air-pollution levels over the surface waters of the bay 
were not higher than over nearby land and why high 
concentrations were observed aloft. We also investigate 
how the Chesapeake Bay breeze affected surface air 
quality and pollutant transport and whether the high 
concentrations of pollutants aloft reached the free tro- 
posphere, as this would have important implications on 
air quality far downwind. 

2. Model description 

In this study, we used the U.S. Environmental Pro- 
tection Agency (EPA) Community Multiscale Air 
Quality (CMAQ; Byun and Schere 2006) model, version 
5.0, which is driven offline by output from the Weather 
Research and Forecasting (WRF; Skamarock et al. 
2008) model, version 3.3, to simulate the state of the 
atmosphere covering the entire months of June and July 
of 2011. The WRF model simulations were performed at 
high resolution (horizontal resolution up to 1.3 km) to 
capture local-scale bay-breeze circulations. The simu- 
lations begin on 24 May 2011 to account for model 
spinup time. 

Meteorological fields are passed to CMAQ at 1-h time 
intervals. Passing the meteorological conditions into an 
air-quality model at a higher temporal resolution or 
running the chemistry online within a meteorological 
model is preferable but requires significantly more 
computational resources (Grell et al. 2004). In our case, 
the demand for computer disk space for passing the 
meteorological fields into the air-quality model at higher 
resolutions or the computational power that would be 
required for running an online meteorological and air- 
quality model was too high to run the models up to 
1.3-km horizontal resolution. Even though performing 
offline high-spatial-resolution air-quality model simu- 
lations with coarse temporal meteorological inputs adds 
some uncertainty to the results by underestimating the 
variability of vertical transport (Grell et al. 2004), this 
approach has been used in a number of previous studies 
to allow for key insights to be gained on how complex 


environments (i.e., topography, coastlines, urban cen- 
ters, and industrial complexes) and/or small-scale pro- 
cesses (i.e., sulfate formation and boundary layer 
venting through fair-weather cumulus clouds and local- 
scale circulations, such as sea breezes) affect air quality 
(Gonsalves et al. 2009; Im et al. 2010; Jimenez et al. 2005, 
2006, 2007; Lee et al. 2008; Loughner et al. 2011; Ortega 
et al. 2009; Parra et al. 2006; Perez et al. 2006; San Jose 
et al. 2007; Stein et al. 2007; Yu et al. 2008). In this paper, 
we use these high-resolution simulations to capture the 
complex geometry of the Chesapeake Bay and accu- 
rately simulate the bay breeze. Previous research has 
shown that coarse-resolution (>5km) model simula- 
tions may not be able to capture this local-scale circu- 
lation (Loughner et al. 2011). 

In addition, recent work has shown that high- 
resolution CMAQ simulations improve the representa- 
tion of vertical transport and boundary layer venting in 
the model (Loughner et al. 2011). It has been suggested 
that CMAQ simulations with resolutions of 12 km un- 
derestimate boundary layer venting on the basis of a pos- 
itive model bias in sulfur concentrations throughout the 
eastern United States (Mueller et al. 2006). Increasing 
convective venting causes pollutants to have longer life- 
times and to be transported greater distances, exacerbat- 
ing air pollution downwind (Cooper et al. 2010). Previous 
studies showed that CMAQ run with a horizontal reso- 
lution of 12 km underestimates interstate transport of 
pollutants (Gilliland et al. 2008; Godowitch et al. 2010). 
Improved meteorological data assimilation in WRF has 
been used to improve horizontal transport (Gilliam et al. 
2012). A more recent study has shown that, as the reso- 
lution increases, vertical mixing increases within the 
boundary layer and between the PBL and the free tro- 
posphere (Loughner et al. 2011). In addition, un- 
certainties and errors in the chemistry and chemical 
lifetime within CMAQ may result in errors in regional 
transport (Gilliland et al. 2008). This study was done by 
comparing CMAQ simulations run with horizontal reso- 
lutions of 13.5, 4.5, 1.5, and 0.5 km (Loughner et al. 2011). 
Results discussed in section 4 below show that the high- 
resolution model simulations are able to capture the 
vertical transport related to the local-scale bay-breeze 
circulation. 

a. Meteorological model 

The WRF model, used to simulate the meteorological 
conditions, is run with one-way nested domains at 36-, 
12-, 4-, and 1.33-km horizontal resolution with di- 
mensions of 149 X 129, 175 X 175, 172 X 220, and 274 X 
307 grid cells, respectively (see Fig. 1 for the model 
domains). All domains have 34 vertical levels from the 
surface to 100 hPa, with 16 levels within the lowest 2 km 
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- 135.0 - 110.0 - 85.0 - 60.0 - 81.0 - 71.0 

Fig. 1. Diagram showing (left) the four modeling domains and (right) the two innermost domains. The four domains 


have horizontal resolutions of 36, 12, 4, and 1.33 km. 


so as to simulate boundary layer processes accurately. 
The domain with 1.33-km horizontal resolution covers 
Washington, Baltimore, and the Chesapeake Bay and is 
at a high enough resolution to simulate local-scale bay- 
breeze events. Model-simulated fields were output 
hourly. The North American Regional Reanalysis is 
used for the model initial and outermost lateral 
boundary conditions. The water surface temperature in 
the model is from the Multiscale Ultrahigh Resolution 
sea surface temperature analysis, which has a resolution 
of 0.01°, or ~1 km. The model is run with the Pleim-Xiu 
surface layer scheme (Pleim 2006) and the Pleim-Xiu 
land surface model (Xiu and Pleim 2001) to calculate 
surface fluxes of heat, moisture, and momentum; the 
Asymmetric Convective Model 2 (ACM2; Pleim 2007) 
to compute mixing in the PBL; the WRF single-moment 
six-class microphysics scheme (WSM-6) to calculate 
water vapor, cloud water, rain, cloud ice, snow, and 
graupel (Hong and Lim 2006); and the Kain-Fritsch 
convective parameterization to simulate subgrid clouds 
and vertical transport within them (Kain 2004). The 
convective parameterization is only turned on in the 
domains with 36-, 12-, and 4-km horizontal resolution. 
Nudging of winds, temperature, and moisture was 
performed for the outermost domain following rec- 
ommendations described in Gilliam and Pleim (2010). 
Gravity waves and vertical velocity damping at the top 
of the model domain are used as described in previ- 
ous work (Loughner et al. 2011; Klemp et al. 2008; 
Skamarock et al. 2008). 

b. Air-quality model description 

The CMAQ model is used to investigate how the bay 
breeze influences pollutant transport and the formation 
of ozone. Here, the model is run with the following user 


options: 1) the Carbon Bond-05 (CB05) gas-phase 
chemical mechanism (Yarwood et al. 2005), 2) CMAQ’s 
fifth-generation modal aerosol model (aero5), 3) ACM2 
for calculating vertical diffusion, 4) the Asymmetric 
Convective Model (ACM; Pleim and Chang 1992) for 
computing convective mixing and containing the het- 
erogeneous chemistry scheme in CMAQ, and 5) the 
Models-3 Dry Deposition (M3DRY) scheme for calcu- 
lating dry deposition (Pleim et al. 2001). Chemical initial 
and boundary conditions come from a Model for Ozone 
and Related Chemical Tracers, version 4 (MOZART-4), 
simulation (Emmons et al. 2010). CMAQ model output is 
saved as hourly averages. The emissions used in the 
CMAQ simulation are described in the appendix. 

3. Observations 

Surface meteorological and trace-gas surface obser- 
vations are used to assess the impact of the bay breeze 
on surface air quality. Both 2-m temperature and 10-m 
wind speeds were observed by the National Weather 
Service, and surface ozone concentrations were mea- 
sured by the Maryland Department of the Environment, 
District of Columbia Department of the Environment, 
Virginia Department of Environmental Quality, and 
Delaware Department of Natural Resources and the 
Environment. Ship observations were made on board 
the National Oceanic and Atmospheric Administration 
(NOAA) Small Research Vessel R8501 (SRVx; NOAA 
Marine Sanctuaries Program) as part of the GeoCAPE- 
CBODAQ field campaign. The ship was equipped 
with a Thermo Environmental Instruments, Inc., model 
49 UV photometric ozone (0 3 ) analyzer. The model 
49 0 3 analyzer determines ambient concentration by 
measuring the attenuation of UV radiation at 254 nm. 
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One -minute-average 0 3 measurements have a 2-sigma 
(std dev) error of 0.8 ppbv. In addition, an NO v [reactive, 
odd nitrogen: nitric oxide (NO), nitrogen dioxide (N0 2 ) 
(NO + N0 2 = NO*), and the compounds produced 
from oxidation of NO. v , such as nitric acid] detector, 
described in Luke et al. (1992), was used on the ship. The 
NOy detector was calibrated using a standard reference 
material (SRM) of n-proly nitrate in the laboratory as 
well as an N0 2 SRM during an in situ calibration. The 
uncertainty of NO v conversion to NO using a molybdenum 
converter heated to 350°C, as performed here, is on the 
order of 15% (Fehsenfeld et al. 1987). 

In situ aircraft observations of 0 3 , NO v , and carbon 
monoxide (CO) made on board the NASA P3-B are 
used to interpret the three-dimensional structure of the 
air pollutants in the atmosphere. Ozone, along with 
reactive nitrogen, is measured with a four-channel 
chemiluminescence instrument with one channel for 
each of NO, N0 2 , NO v , and 0 3 . The NO is measured by 
adding reagent 0 3 to the sample flow stream, producing 
excited N0 2 , which emits photons in front of a dry-ice- 
cooled photomultiplier tube. Photons are counted to 
provide the fundamental signal. N0 2 is measured in 
a separate channel via photolytic conversion to NO 
using 400-nm light-emitting diodes, followed by de- 
tection as NO. In a similar way, NO v is measured via 
catalytic conversion of reactive nitrogen species to NO 
through reaction with CO on a solid gold surface 
heated to 300°C. Ozone is measured using the same 
reaction by the addition of reagent NO to the sample 
flow. The 1-s detection limit for NO is lOpptv. It is 15 
pptv for NO? and NO v and 0.1 ppbv for 0 3 . At high 
mixing ratios, the overall uncertainty for NO and N0 2 is 
10%, and it is 5% for 0 3 . 

The CO, along with methane (CH 4 ) and nitrous oxide 
(N 2 0), is measured in situ using a three-channel tun- 
able diode laser absorption instrument known as 
“DACOM” (Sachse 1987, 1988). Dried ambient air is 
drawn continuously through a low-pressure multipass 
optical cell where near-infrared diode laser beams 
measure concentrations using line-locked wavelength 
modulation spectroscopy. Accuracy is maintained through 
frequent introduction of calibration gases whose con- 
centrations are determined separately by the appropri- 
ate groups at the NOAA Earth Systems Research 
Laboratory/Global Monitoring Division. Data are re- 
ported at one sample per second. The CO observations 
have a precision of 1 ppbv and an uncertainty of 1.5%. 

4. Results 

Here, we discuss how a bay breeze affected surface air 
quality and pollutant transport within and between the 


PBL and the free troposphere on 11 July 2011 using 
observations and model results from the 1.33-km- 
horizontal-resolution domain. Prior to the start of the 
bay breeze, overnight and early-morning winds were 
from the south, transporting pollutants over Washington 
and Baltimore, northward into Pennsylvania (Fig. 2). The 
southerly winds prevented surface air pollution over the 
Washington-Baltimore urban corridor from being ad- 
verted over the Chesapeake Bay. Both observations and 
model calculations showed that the bay-breeze event 
began early on 11 July along the western coastline of the 
Chesapeake Bay. By 1100 eastern daylight time (EDT), 
a strong virtual potential temperature gradient along the 
western coastline of the bay initiated this change in wind 
direction and the bay breeze began to form (Fig. 2). We 
point out the virtual potential temperature gradient, 
which is a proxy for the density gradient that drives the 
bay-breeze circulation. The bay breeze persisted and 
grew in strength throughout the day. The observations 
and WRF model simulation show that the bay-breeze 
convergence zone penetrated inland into the urban cor- 
ridor, with strong surface convergence present by late 
afternoon (Fig. 2). A vertical cross section of virtual po- 
tential temperature and wind velocities at 1400 EDT 
shows winds converging and upward vertical motion at 
the bay-breeze convergence zone near Padonia, Mary- 
land (Fig. 3). 

The model simulation for 11 July is compared with 
ground-, aircraft-, and ship-based observations. Mod- 
eled and observed ozone are compared by calculating 
mean bias, normalized mean bias, root-mean-square 
error, and normalized mean error, as defined by Eder 
and Yu (2006) for maximum 8-h average at ground- 
based monitoring sites, 15-s-average data observed by 
the NASA P-3B below 5 km AGL, and hourly-average 
data measured on board the ship (Table 1). The ozone- 
model bias can partially be explained by uncertainties in 
the emissions and chemistry in the model. The model 
overestimated N O v concentrations, suggesting emis- 
sions are too high (Anderson et al. 2014, manuscript 
submitted to Atmos. Environ.). Modeled NO v concen- 
trations were 25% larger than ship-observed concen- 
trations and 288% larger than NASA P-3B aircraft 
observations on 11 July. The CMAQ model over- 
estimated NO y concentrations in comparison with ship 
observations during most of the days of the CBODAQ 
field campaign (Goldberg et al. 2014). The over- 
estimation of NO v over other U.S. cities has been seen 
during other modeling studies as well (Brioude et al. 
2013; Yu et al. 2012). The positive bias in NO v is also 
explained by uncertainties in conversion rates of NO, 
(i.e., NO v — N O r ) to N0 2 within the CB05 chemical 
mechanism used in CMAQ as described by Goldberg 
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Fig. 2. (left) Observed 2-m temperature and 10-m wind velocity, (center) WRF-diagnosed 2-m temperature and 10-m wind velocity, and 
(right) WRF-simulated virtual potential temperature in the first model layer (about 11 m AGL) and 10-m wind velocity on 11 Jul 2011 at 
(top) 0700, (middle) 1100, and (bottom) 1400 EDT. The WRF results are from the 1.33-km-horizontal-resolution domain. One of every 14 
grid cells has a wind vector plotted. The bay-breeze circulation along the northern half of the western coastline of the bay began to form at 
1100 EDT when winds over the water changed from calm or southwesterly to south-southeasterly while winds over land remained 
southerly or south-southwesterly (see middle row). The bay-breeze convergence zone pushed farther inland by 1400 EDT (see bottom 
row). The letters W and B in the top-left panel show the locations of Washington and Baltimore, respectively, and the letter P in the 
middle-left panel denotes the location of Padonia. The blue lines in the middle-right and bottom-right panels show the location of the bay- 
breeze front. The bay breeze was identified by noting the change in surface wind direction associated with the large virtual potential 
temperature gradient, surface convergence, and vertical motion at the bay-breeze front, as shown in Fig. 3, below. 


et al. (2014). We focus on the role of the bay breeze in 
transporting air pollutants out of the boundary layer as 
shown in the P-3B observations and the model simula- 
tion (see below). 


Observations from ground-based monitoring sen- 
sors showed localized areas of high values of maximum 
8-h-average surface ozone concentrations near the 
bay-breeze convergence zone, which was within the 
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Fig. 3. Vertical cross section of virtual potential temperature and 
explicit wind velocities at 1400 EDT in the lowest 1 km between 
77.1° and 76.1°W at 39.46°N. One of every six grid cells in the 
horizontal direction has a wind vector plotted. This cross section 
passes through Padonia, which is located at 39.46°N, 76.63°W. 

Washington-Baltimore metropolitan area (Fig. 4). As 
the bay-breeze convergence front penetrated within the 
urban corridor, the surface winds caused urban emis- 
sions on both sides of the bay-breeze front to converge. 
Farther away from the convergence zone, cleaner air 
over the water and over rural and suburban areas pen- 
etrated into urban areas, resulting in a large horizontal 
gradient in ozone concentrations with the highest 
concentrations located near the bay-breeze conver- 
gence zone in the urban corridor. The model positive 
bias in ozone (Fig. 4) is due to uncertainties in the 
emissions and chemistry as discussed above and in 
Goldberg et al. (2014). 

The southerly winds prevented pollutants from 
Washington and Baltimore from being transported over 
the bay during the overnight and early-morning hours 
prior to the development of the bay breeze (Fig. 2), 
preventing higher surface ozone concentrations over the 
Chesapeake Bay than over the adjacent land (Fig. 4). 
The locations of the ship throughout 11 July as part of 
the GeoCAPE-CBODAQ field campaign are shown in 
Fig. 5. The model positive bias exists over water (Fig. 6) 
as well as over land (Fig. 4); both model and observations 
showed lower surface ozone values over the estuarine 
waters than over the land, however. 

The NASA P-3B aircraft performed in situ observa- 
tions over the region and performed spirals sampling the 
air in the vertical direction to obtain a three-dimensional 
view of air-pollution concentrations (Fig. 5). Both the 


Table 1. The 11 Jul 2011 CMAQ 1.33-km-horizontal-resolution- 
domain simulated mean bias (MB), normalized mean bias (NMB), 
root-mean-square error (RMSE), and normalized mean error 
(NME), as defined by Eder and Y u (2006), for maximum 8-h-average 
0 3 at ground-based monitoring sites, 15-s average P-3B 0 3 mea- 
surements within and above the PBL, and hourly-average ship 0 3 
observations. 



MB 

NMB 

(%) 

RMSE 

NME 

(%) 

Ground sites max 8-h avg 

16.3 

29.3 

18.0 

29.3 

P-3B hourly avg within 

9.7 

12.8 

16.6 

16.9 

PBL 

P-3B hourly avg above 

-0.66 

-0.97 

14.6 

18.2 

PBL 

Ship hourly avg 

9.9 

17.6 

10.6 

17.6 


model and aircraft observations reveal that ozone con- 
centrations throughout the PBL were highest near the 
bay-breeze convergence zone near Beltsville (Mary- 
land) and Padonia (Fig. 7). Rapid mixing in the PBL 
caused high air-pollution concentrations near the sur- 
face to be transported throughout the boundary layer. 
Winds at 850 hPa, or about 1.5 km, within the boundary 
layer near the bay-breeze convergence zone and in the 
free troposphere downwind of this region were south- 
westerly (Fig. 8). These upper-level winds transported 
pollutants near the top of the boundary layer near the 
bay-breeze convergence zone toward Aldino, Edge- 
wood, Essex, and Fair Flill (all in Maryland), where 
a lower PBL height was present than in the region of the 
convergence zone (Fig. 7) because of southerly or 
southeasterly surface winds from the bay advecting 
cooler air from the bay over the four observational sites 
(Fig. 2). Pollution near the top of the boundary layer at 
the bay-breeze convergence zone was transported 
downwind where it was horizontally advected out of the 
PBL and into the free troposphere. High ozone con- 
centrations were observed and simulated in the free 
troposphere over Aldino, Edgewood, Essex, and Fair 
Hill (Fig. 7). Note that the bay breeze on this day oc- 
curred on the western side of the bay and not the 
northern end. The northern coastline of the bay experi- 
enced offshore winds before and after the bay breeze 
formed. Areas near the northern coastline of the bay 
(i.e., Aldino and Edgewood) were affected by the bay 
breeze aloft, however. The upper levels over this area 
were downwind of the bay-breeze front, resulting in high 
air-pollution levels aloft (Fig. 7). In a similar way, 
comparisons between the NASA P-3B and CMAQ re- 
veal high CO concentrations at the surface and through- 
out the boundary layer near the bay-breeze convergence 
zone and in the free troposphere downwind of the bay- 
breeze convergence zone (Fig. 9). 
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Fig. 4. (left) Observed and (right) CMAQ-simulated maximum surface 8-h-average ozone on 11 Jul 2011. CMAQ 
results are from the 1.33-km-horizontal-resolution domain. 


A time series of CMAQ-simulated ozone maps at the 
surface reveals that the highest surface concentrations 
in the Washington and Baltimore metropolitan areas 
were located near the bay-breeze convergence zone af- 
ter the early formation of the bay breeze (Fig. 10; 1200- 
1800 EDT). In addition, ozone maps at about 1.5 km 
AGL show that increased ozone concentrations aloft 
originated near the bay-breeze convergence zone during 
the early afternoon (Fig. 10; 1400-1800 EDT). The high 
levels of ozone and its precursors were then transported 
downwind to the east-northeast. It can be seen that 
CMAQ simulated high ozone concentrations aloft over 
areas with low surface ozone concentrations north of the 
Chesapeake Bay (Fig. 10; model simulations for 1600 
and 1800 EDT), which is also shown in CMAQ and P-3B 
comparisons (Fig. 7). 

A back trajectory was calculated to investigate 
where the high pollution concentrations at 1.5 km over 
Edgewood originated, and a forward trajectory was 
calculated to determine where this plume was trans- 
ported (Fig. 11). The domain with 4-km horizontal res- 
olution was used because it has more area than the 
domain with 1.3-km horizontal resolution yet the 
resolution was high enough to capture the local-scale 
bay-breeze circulation. Kinematic three-dimensional 
backward and forward trajectories from Edgewood ini- 
tialized at 1.5 km AGL at 1600 EDT show that the air 
parcels originated in the boundary layer west of the bay- 
breeze front at 1400 EDT, traveled northeast over 
Edgewood at 1600 EDT above the PBL, and entered 
southern New England by 0200 EDT the next day 
(Fig. 11). The trajectory exited the 4-km domain after 


0300 EDT. Cross sections of ozone and carbon monox- 
ide along the trajectories show high pollutant concen- 
trations throughout the boundary layer between 1400 
and 1500 EDT, which was when the back-trajectory path 
was within the boundary layer near the bay-breeze front 
(Figs. 12 and 13). The cross sections show increased 
pollution aloft (from 500 m to 2 km in altitude) from 



Fig. 5. Map showing the location of the ship dock and the six 
stations (1-6) to which it went on 11 Jul 2011. Also shown is the 
flight path of the NASA P-3B on 11 Jul 2011, the locations where 
the plane spiraled over monitoring sites, and the location of 
Washington and Baltimore. The P-3B made two clockwise circuits, 
as depicted by the arrows. 
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Fig. 6. Observed (black) and CMAQ-simulated (red) surface 
ozone concentrations from the ship on the Chesapeake Bay on 11 
Jul 2011. The black lines and the numbers 1-6 depict when the ship 
was at each of the six stations shown in Fig. 5. CMAQ results are 
from the 1.33-km-horizontal-resolution domain. 



1600 to 0300 EDT, which was when the trajectory was 
downwind of the bay-breeze convergence zone in the 
free troposphere (Figs. 12 and 13). Peak 0 3 and CO 
mixing ratios descended in altitude overnight to ~1 km, 
making it likely that this plume mixed down to the sur- 
face as the boundary layer grew in depth on 12 July 


outside the area of the 4-km domain, affecting surface 
air quality in New England. 

5. Discussion 

We demonstrated how a bay breeze caused locally 
high surface ozone concentrations on 11 July 2011 in 
Maryland. The bay-breeze convergence zone located 
within the Washington-Baltimore metropolitan area 
caused high air-pollution concentrations on both sides of 
the bay-breeze front to converge. Cleaner surface air 
over the water and over rural and suburban areas was 
transported toward the bay-breeze convergence zone, 
which resulted in a large horizontal gradient in ozone 
concentrations. The highest observed maximum 8-h- 
average ozone reached 75 ppbv at the bay-breeze con- 
vergence zone in Beltsville. 

A key ingredient in exacerbating air-pollution levels 
further that was missing in this case was the recirculation 
of air pollutants as described by Loughner et al. (2011) 
and summarized here. When synoptic-scale surface 
winds are westerly, the wind direction over the water 
must change almost 180° for a bay breeze to impact the 
western shore of the Chesapeake Bay. Offshore winds 
transport pollutants over land across the coastline out 
over the bay waters. As the winds reverse directions 
over the water, stagnant conditions ensue, resulting in 


5 


4 



0 40 45 50 55 60 65 70 75 80 85 90 

Ozone (ppbv) 


Fig. 7. CMAQ simulated (background) and observed (overlay) ozone concentrations along a flight track on 11 Jul 
2011. The white line shows the location of the top of the boundary layer as calculated by the WRF model. The black 
letters at the bottom of the figure — Be, Pa, Fa, Al, Ed, Es, and CB — stand for the spiral locations over Beltsville, 
Padonia, Fair Hill, Aldino, Edgewood, Essex, and the Chesapeake Bay, respectively. CMAQ results are from the 
1.33-km-horizontal-resolution domain. 




1706 


JOURNAL OF APPLIED METEOROLOGY AND CLIMATOLOGY 


Volume 53 



288 289 290 291 292 293 294 295 296 10 m/s 

Temperature (K) 10 m wind speed 


Fig. 8. WRF temperature and wind velocities at 850 hPa, which is at an altitude of about 1.5 km, at (left) 1400 and 
(right) 1600 EDT. WRF results are from the 1.33-km-horizontal-resolution domain. One of every 14 grid cells has 
a wind vector plotted. The red letters in the left panel — Be, Pa, Fa, Al, Ed, and Es — stand for the spiral locations over 
Beltsville, Padonia, Fair Hill, Aldino, Edgewood, and Essex, respectively. 


pollutants accumulating over the water. After the bay 
breeze forms, the local-scale winds transport the pollut- 
ants over the bay back inland to the convergence zone 
where they converge with freshly emitted pollutants 
over land. On 11 July 2011, however, the synoptic-scale 
winds were parallel to the western coastline of the bay. 


Southerly winds during the overnight and early-morning 
hours prevented surface air pollution in the Washington 
and Baltimore metropolitan area from being transported 
across the coastline prior to the onset of the bay breeze. 

The southerly winds on 11 July played a role in the 
formation and strength of the bay breeze. The onset of 
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Fig. 9. As in Fig. 7, but for carbon monoxide. 
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Fig. 10. CMAQ-simulated ozone concentrations (left) at the surface and (right) at —1.57 km 
AGL at (top) 1200, (top middle) 1400, (bottom middle) 1600, and (bottom) 1800 EDT from the 
1 .33-km-horizontal-resolution domain. 
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Fig. 11. Horizontal path of a back trajectory from 1.5 km 
AGL over Edgewood (black) from 1600 to 1400 EDT (red) and 
a forward trajectory from 1.5 km AGL over Edgewood from 
1600 to 0200 EDT (blue) using the 4-km-horizontal-resolution 
WRF simulation. 

a bay-breeze circulation begins earlier, when the 
synoptic-scale winds are parallel to the coastline with 
land to the left (lower pressure over land), which is 
called a “corkscrew” bay breeze, as described in detail 
by Steele et al. (2013) and summarized here. This was 
the case on 11 July along the western shore of the 
Chesapeake Bay. A corkscrew bay breeze can be 


pictured as winds following the coastline in a corkscrew 
pattern with rising air over land and sinking air over 
water. Prior to the corkscrew bay breezes forming, 
stronger surface friction over land than over water re- 
sults in divergence and descending motion at the 
coastline. This descending motion helps to initiate the 
bay-breeze circulation, allowing the local-scale circulation 
to form earlier, with a weaker surface virtual potential 
temperature gradient, and under stronger early-morning 
winds than if the synoptic-scale winds were offshore. 
Likewise, if the synoptic-scale winds were parallel to the 
coastline with the land to the right (higher pressure over 
land), as on the eastern shore of the Chesapeake Bay on 
11 July, differences in surface friction in the early 
morning would cause convergence and rising motion at 
the coastline, which would make it more difficult for 
a bay breeze to develop (Steele et al. 2013). A bay 
breeze that forms under these conditions is called 
a backdoor bay breeze, and it would be weaker and form 
later in the day than corkscrew bay breezes under similar 
wind speeds and virtual potential temperature gradients 
(Steele et al. 2013). As opposed to other bay-breeze 
events, corkscrew bay-breeze events do not magnify the 
trace-gas mixing ratios in air-pollution events through 
recirculation of surface air pollution, but they do amplify 
the duration of the air-pollution events through pro- 
longed strong surface convergence. 

In addition, we see that strong, prolonged bay-breeze 
events can transport a large amount of air pollution out 
of the PBL and into the free troposphere. The lifting 
motion at the sea- and bay-breeze convergence zones 
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Fig. 12. Vertical cross section of CMAQ-simulated ozone concentrations along the trajec- 
tories in Fig. 11 from the 4-km-horizontal-resolution domain. The white line shows the location 
of the top of the boundary layer as calculated by the WRF model from the 4-km-horizontal- 
resolution domain. 
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Fig. 13. As in Fig. 12, but for carbon monoxide. 


has been shown to initiate thunderstorms and severe 
weather (Kingsmill 1995), as well as to transport pol- 
lutants aloft (Loughner et al. 2011). Little research has 
been performed to show how bay-breeze circulations 
affect boundary layer venting. This process is significant 
in that pollutants transported from the boundary layer 
to the free troposphere affect air quality farther down- 
wind. Once in the free troposphere, pollutants have lon- 
ger lifetimes and are susceptible to long-range transport. 
These pollutants in the free troposphere can eventually 
be transported back into the PBL and subside to the 
surface, affecting air quality, human health, and ecosys- 
tem processes far away from the pollution source 
(Cooper et al. 2010). In this study, we show that the 
pollution layer remained aloft as it traveled into southern 
New England. We did not determine where this layer 
eventually reentered the boundary layer and mixed down 
to the surface to affect air and water quality because it 
was transported out of our domain. 

6. Conclusions 

In this study, WRF and CMAQ model simulations are 
analyzed alongside ground-based, aircraft, and ship- 
board observations to investigate the role of the Ches- 
apeake Bay breeze on surface air quality, pollutant 
transport, and boundary layer venting. Southerly winds 
during the overnight and early-morning hours prevented 
pollutants from the Washington and Baltimore area 
from being transported out over the surface waters of 
the Chesapeake Bay. Therefore, the pollutants were not 
recirculated over the water and then back inland once 
the bay breeze formed to converge with freshly emitted 
pollutants over land near the bay-breeze convergence 


zone. Southerly synoptic-scale winds are, however, fa- 
vorable for producing a strong corkscrew bay breeze 
that forms early in the day along the western coastline of 
the Chesapeake Bay. Results show that a strong and 
prolonged bay breeze caused locally high surface ozone 
concentrations near the bay-breeze convergence zone, 
which penetrated into the urban corridor. The onset of 
the bay breeze occurred early, by 1100 EDT, resulting in 
an extended period of convergence of air pollution at 
the bay-breeze front. 

This strong bay-breeze circulation also resulted in 
a significant amount of air pollution being transported 
out of the PBL and into the free troposphere, which has 
implications for regional air quality. Pollutants trans- 
ported from the PBL to the free troposphere gain a longer 
lifetime and are susceptible to long-range transport. These 
pollutants can then subside back into the PBL, having an 
impact on surface air quality, human health, and ecosys- 
tem processes far away from their emissions sources. 
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APPENDIX 

Emissions 

Anthropogenic emissions input files for the CMAQ 
model are created with the Sparse Matrix Operator 
Kernel Emissions (SMOKE) modeling system (Houyoux 
and Vukovich 1999). Because a 2011 emissions inventory 
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is not yet available, the closest projected inventory was 
used. We used a projected 2012 emissions inventory, 
which is grown from the 2005 National Emissions In- 
ventory (NEI) to include estimated emissions changes 
due to growth and emissions controls that were to be 
implemented by 2012 (EPA 2011). Annual projected 
point and countywide area emissions are temporally 
distributed on the basis of time of day, day of week, and 
season using temporal surrogates from the EPA. Mobile 
emissions estimates from cars, trucks, and motorcycles 
are computed with the Motor Vehicle Emission Simu- 
lator (MOVES; EPA 2012). Point sources are vertically 
distributed on the basis of the meteorological condi- 
tions, stack height, and the temperature and velocity of 
the emissions exiting the stack. Countywide area emis- 
sions are horizontally distributed on the basis of land use 
from spatial surrogates provided by the EPA for the do- 
mains with 36- and 12-km horizontal resolution. Spatial 
surrogates for the 4- and 1.33-km domains are generated 
from the Multimedia Integrated Modeling System 
(MIMS) Spatial Allocator (Eyth and Brunk 2005) using 
data from the 2000 U.S Census, National Land Cover 
Characteristics Data, and other spatial sources available 
from EPA’s Emissions Modeling Clearinghouse. 

Biomass burning emissions are obtained from the 
Fire Inventory from NCAR, version 1.0 (FINNvl; 
Wiedinmyer et al. 2011). This inventory estimates 
biomass-burning emissions with a resolution of 1 km on 
the basis of fire hot spots, area burned, land-cover maps, 
and biomass consumption estimates from Moderate 
Resolution Imaging Spectroradiometer (MODIS) sat- 
ellite products (Wiedinmyer et al. 2011). In the ideal 
case, the emissions would be dispersed in the vertical 
direction on the basis of the meteorological conditions 
and the intensity of the fire. A plume-rise algorithm for 
vertical dispersion of fire emissions in CMAQ has not 
been developed, however, and therefore the emissions 
are placed in the lowest level of the model. Plumes of 
biomass-burning emissions were not observed in the 
Washington and Baltimore metropolitan area during 
the time of interest to this study, and therefore placing 
the fire emissions near the surface did not have any 
impact on the model results presented here. Airborne 
observations made on board the NASA P-3B during the 
DISCOVER-AQ field campaign showed negligible 
wildfire impact in the region. Measurements of aceto- 
nitrile, a tracer for biomass burning, revealed that 
no smoke plumes were present in the Washington- 
Baltimore metropolitan region during July. Observed 
acetonitrile observations on 11 July were low, with 
minimum, median, and maximum concentrations of 
0.102, 0.144, and 0.198 ppbv, respectively. Observed 
acetonitrile concentrations over the western United 


States from long-range transported biomass-burning 
plumes originating from Asia have been observed to 
be —0.5 ppbv (de Gouw et al. 2004). 

Biogenic and lightning emissions are calculated online 
within the CMAQ model. The Biogenic Emissions In- 
ventory System (BEIS) calculates biogenic emissions on 
the basis of the meteorological output from the WRF 
simulation and land-use data describing the vegetation 
(Vukovich and Pierce 2002). Lightning NO* emissions 
are calculated from lightning-flash count data from the 
National Lightning Detection Network and convective 
precipitation as calculated in WRF (Allen et al. 2012). 
Convective precipitation is calculated within WRF’s 
convective parameterization, and, because no convective 
parameterization is turned on in the finest horizontal 
domain (1.33-km resolution), the lightning emissions 
from the 4-km domain are used in the 1.33-km domain. 
This method may result in the placement of lightning 
emissions in areas where deep convection is not located 
in the 1.33-km domain since parameterized convection 
in the 4-km domain and resolved convection in the 

1.33-km domain could behave differently. Misplaced 
lightning emissions do not affect the results during the 
11 July bay-breeze scenario discussed in this paper, 
however. No lightning emissions are placed in the 

1.33-km domain on 10 and 11 July. The area of the 

1.33-km domain is small, and so the effects of lightning 
emissions placed in the 1.33-km domain on 9 July were 
transported well out of the area of interest. Therefore, the 
only lightning emissions affecting the 1.33-km domain, if 
any, during the 11 July case study presented in our paper 
are emissions calculated online in the outer domains that 
are transported into the innermost domain. 
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